Ab initio Derivation of Low-Energy Model for k-ET Type Organic Conductors 
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We derive effective Hubbard-type Hamiltonians of k-(ET)2X, using an ab initio downfolding technique, for 
the first time for organic conductors. They contain dispersions of the highest occupied Wannier-type molecular 
orbitals with the nearest neighbor transfer i~0.067 eV for a metal X=Cu(NCS)2 and 0.055 eV for a Mott 
insulator X=Qi2(CN)3, as well as screened Coulomb interactions. It shows unexpected differences from the 
conventional extended Hiickel results, especially much stronger onsite interaction U ~0.8 eV (U /t~ 12-15) than 
the Hiickel estimates (?7/t~7-8) as well as an appreciable longer-ranged interaction. Reexamination on physics 
of this family of materials is required from this realistic basis. 
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Organic conductors with ET molecules, (ET) 2 X with a 
number of choices of anions X, offer a variety of prototypical 
behaviors of strongly correlated electron systems with two- 
dimensional (2D) anisotropics [1]. Examples range from cor- 
related metals with superconductivity at low temperatures to 
Mott insulators either with a quantum spin liquid or with anti- 
ferromagnetic, charge-density or spin-Peierls orders. Intrigu- 
ing Mott transitions are also found. They are all in front of 
recent active research for unconventional quantum phases and 
quantum critical phenomena in nature, while their essences of 
physics are still under strong debates. 

In particular, an unconventional nonmagnetic Mott- 
insulating phase is found near the Mott transition in the K-type 
structure of ET molecules, X=Cu2(CN)3 referred to as k-CN, 
where no magnetic order is identified down to the tempera- 
ture T=0.03 K, four orders of magnitude lower than the an- 
tiferromagnetic spin-exchange interaction J~250 K yfl. The 
emergence of the quantum spin liquid near the Mott transi- 
tion has been predicted in earlier numerical studies ill, 
while the full understanding of the spin liquid needs more 
thorough studies. It is also crucially important to elucidate 
the real relevance of the theoretical findings to the real k-ET 
compounds. Most of numerical [6] and theoretical [0] stud- 
ies have also been performed for a simplified single-band 2D 
Hubbard model based on an empirical estimate of parameters 
combined with extended Hiickel calculations ISllSD- A more 
realistic description of k-ET compounds is certainly needed 
beyond the empirical model. 

Another fundamental finding achieved in this series of 
compound is the unconventional Mott transition found for 
X=Cu[N(CN)2]Cl under pressure [10]. The novel univer- 
sality class of the Mott transition is in good agreement with 
the marginal quantum criticality at the meeting point of the 
symmetry breaking and topological change [IT, 12, Til [Till . 



Because of its significance to the basic understanding on the 
physics of quantum criticality, the relevance of theoretical 
concept to the experimental observation needs to be further 
examined on the realistic and first-principles grounds. Fur- 
thermore, an unconventional superconductivity is found in 



the metallic sides of these compounds at low temperatures 
(T<T C ~10-13K), where the mechanism is not clearly under- 
stood yet 015[|l6ll . These outstanding properties of k-ET com- 
pounds encourage systematic studies based on realistic basis. 
As mentioned above, however, the first-principle studies are 
limited IU7I1 and most of the studies so far were performed 
using the empirical models inferred from the Hiickel studies. 

The purpose of this letter is to derive ab initio effective 
Hamiltonian of the real k-ET compounds to serve in estab- 
lishing a firm basis for studies on the unconventional phe- 
nomena and open issues. The present work is the first chal- 
lenge aiming at ab initio model construction of organic con- 
ductors containing a large number of atoms with four com- 
plex ET molecules in a unit cell. We derive models of two 
contrasting compounds, spin-liquid k-CN and superconduct- 
ing compound Jf=Cu(NCS) 2 abbreviated as k-NCS (TH], to 
get insights into the whole series of K-(EThX compounds 
from metals to Mott insulators. Our ab initio results indicate 
a substantial difference from the previous simple extended 
Hiickel results [HISS; in particular, a large ratio of the onsite 
Coulomb repulsion to the nearest neighbor transfer f7/t~12 
even in metallic k-NCS compared to ~7 in the Hiickel esti- 
mate requires reexamination of the model studies. 

A reliable first-principle framework has been proposed for 
the methods to derive effective low-energy Hamiltonians of 
real materials. This treatment consists of ab initio density- 
functional calculations of the global electronic band struc- 
ture and a subsequent downfolding procedure by elimina- 
tion of degree s of freedom far away from the Fermi level 
1 19 , 2(1 21 , 22 1 . Recently, the technique was combined with a 
framework of maximally localized Wannier orbital (MLWO) 
lEl H. The MLWO Q has a distinct computational ad- 
vantage, because it enables construction of localized orbitals 
even for the present complex molecular solid 02611 . Below 
we will demonstrate how the Wannier-framework downfold- 
ing and the derivation of effective models have successfully 
been worked out. 

Now, we consider ab initio derivations of a single-band ex- 
tended Hubbard Hamiltonian describing the low-energy elec- 
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tronic property of /i-(EThX. For the system, the basis of 
the Hamiltonian is the Wannier functions associated with an- 
tibonding states of the highest occupied molecular orbitals 
(HOMOs) of ET molecules forming a dimer. The explicit 
form of this Hamiltonian is given by 



(1) 



up l] 



where a\ a (aw) is a creation (annihilation) operator of an elec- 
tron with spin a in the Wannier orbital localized at the ith ET 
dimer. The parameters are written by 



U 



Wh\<l>3 



(2) 



with \<pi) — a||0) and Ho being the one-body part of TL The 
Vij values are screened Coulomb integrals in the Wannier or- 
bital, expressed as 



drdr'#(r)&(r)W(ry)#(r')&(r'), (3) 



where W(r, r') is a screened Coulomb interaction. The Vy 
at i = j corresponds to onsite Hubbard repulsion U. Here, W 
is practically calculated in the reciprocal space as 
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where q is a wave vector in the first Brillouin zone and G 
is a reciprocal lattice vector. The crystal volume is fl and 
£ G ' G ,(q) is the inverse dielectric matrix which is related to 
the irreducible polarizability \ by £gg' (<z) = <5gg' — v(q + 
G)xGG'{q)> where v(q)=4ir /£l\q\ 2 is the bare Coulomb in- 
teraction. The polarization matrix xgc (?) with a random 
phase approximation (RPA) is calculated as 

XGG'(Q)-E E^ +q |e J ( q+G >|^ k ) 

k a/3 

jak+q - //3k 



x(V>/3k|f 



i(q+G)r 
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where ijj a k, £ a k, and / Q k are the Bloch state, its energy, and 
the occupancy, respectively. The calculation of x is straight- 
forward with Eq. (O but there is a notice when we consider 
the model construction. The screened interaction put in the 
Hubbard model [Vij in Eq. ([T]i] should not include screening 
formed in a target band of the model lfl9i I20I1 . This screening 
should be considered at the step of solving the effective model 
and, at the downfolding stage, we must exclude the target- 
band screening effects in the calculation of Vij . In RPA, this 
exclusion is easily implemented, since the RPA polarization 
in Eq. <(5j is written in terms of the sum of the band pairs a/3 
associated with individual transitions. We thus generate the 
polarization function with cutting the transitions between the 
target bands and then evaluate the screened interaction W by 
using this polarization function. Finally, we compute the Vij 



parameters as the Wannier matrix elements of the W inter- 
action. The practical details for the matrix evaluation can be 
found in Ref . |24j. 

Our ab initio calculations were performed with Tokyo Ab 
initio Program Package B27I1 . With this program, density- 
functional calculations H28I1 with the generalized-gradient- 
approximation (GGA) exchange-correlation functional ll29ll 
were performed using a plane-wave basis set and the 
Troullier-Martins norm-conserving pseudopotentials I30I1 in 
the Kleinman-Bylander representation 03111 . The energy cut- 
off in the band calculation was set to 36 Ry and a 3x5x5 
fc-point sampling was employed. The experimental crystal- 
structure data were taken from Ref. 18 for k-NCS and Ref.[l6l 
for k-CN. The polarization function was expanded in plane 
waves with an energy cutoff of 3.5 Ry and the total number 
of bands considered in the polarization calculation was set to 
634 for k-NCS and 612 for k-CN. [The total number of the 
occupied (unoccupied) states is 234 (400) for k-NCS and 242 
(370) for k-CN.] The Brillouin-zone integral over wavevec- 
tor was evaluated by the tetrahedron method. The additional 
terms in the long- wavelength polarization function due to non- 
local terms in the pseudopotentials was explicitly considered 
following Ref. [32l The singularity of the Coulomb interaction 
in the q^Q limit, in the evaluation of the elements Vij, was 
carefully handled on the manner described in Ref. |32l 

We show in Fig.[T|aZ? initio GGA band structures of k-NCS 
(a) and k-CN (b). The ab initio results clearly show that the 
band near the Fermi level (energy zero) is totally isolated from 
other bands, thus justifying employing this band as the target 
band of the extended Hubbard model. The entangled bands 
below —0.5 eV and above +1 eV are associated with deeper 
and upper states of the ET molecules, as well as the anion elec- 
tronic structures. We find that the bandwidth of k-NCS esti- 
mated as 0.56 eV is larger than that of k-CN (0.45 eV), which 
generates a discernible difference in transfer integrals of the 
two compounds (see below). We also note that the ab initio 
dispersion of the target band is basically similar to the result 
of the extended Hiickel model but our bandwidth is somewhat 
smaller than the Hiickel results of 0.57 eV for k-NCS and 0.50 
eV for k-CN lH. 



TABLE I: List of the parameters in the single-band extended Hub- 
bard Hamiltonian in Eq. (0 for k-(ET)2X. 



X=Cu(NCS) 2 X=Cu 2 (CN) 3 



t a (meV) 


-64.8 


-54.5 


t b (meV) 


-69.3 


-54.7 


t e (meV) 


44.1 


44.1 


t d (meV) 


-11.5 


- 6.8 


U(eV) 


0.83 


0.85 


X (eV^ir 1 ) 


5.1 


5.0 



Figure Q](c) visualizes our calculated maximally localized 
Wannier orbitals for k-NCS, for the isolated band discussed 
above. An initial guess for generating the Wannier function is 
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FIG. 1: Calculated ab initio GGA band structures (red line) of 
k-(BEDT-TTF) 2 Cu(NCS) 2 (a) and k>(BEDT-TTF) 2 Cu 2 (CN) 3 (b). 
The zero of energy is the Fermi level. The blue dotted disper- 
sions are obtained by the four transfer parameters listed in TABLE I. 
(c) Calculated maximally localized Wannier functions of k-(BEDT- 
TTF) 2 Cu(NCS) 2 . The amplitudes of the contour surface are +1.5/^/v 
(blue) and —\.5ly/v (red), where v is the volume of the primitive 
cell. S, C, H, N, and Cu nuclei are illustrated by green, yellow, sil- 
ver, blue, and red spheres, respectively, (d) Schematic diagram for 
transfer network in the triangular lattice. 



the p-type Gaussian sitting on the dimer center. The resulting 
plot for the Wannier orbital clearly exhibits that the basis state 
of the low-energy effective model of the k-NCS compound is 
the antibonding states consisting of a linear combination of 
the ET HOMOs in the dimer [33]. We note that there is no 
discernible difference between the MLWO of k-NCS and that 
of k-CN. 

The transfer integrals in Eq. (|2]i calculated as matrix ele- 
ments of the Kohn-Sham Hamiltonian in the Wannier orbital 
are listed in TABLE I. We list only the four parameters (t a , tt, 
t c , and td) and the absolute value of other transfers are all less 
than 5 meV. The schematic diagram for the transfer network 
is illustrated in Fig. Q] (d). The band dispersion calculated 
with the four transfers is given as blue dots in Figs. Q] (a) and 
(b) and, from the figures, we see that the original band struc- 
ture is well reproduced with these four transfers. Consider- 
ing t a ^ti,=t (nearest-neighbor transfer), t c =t' (next-neighbor 
transfer), and id~0, a rough estimate for geometrical frustra- 
tion |t7t| is 0.66 for k-NCS and 0.80 for k-CN, thus indicat- 
ing that the degree of frustration in k-CN is larger than that in 
k-NCS, while the estimates are still far from the perfect situa- 
tion |i'/i|=l inferred by the extended Hiickel result for k-CN 



dielectric function for k-NCS (a) and k-CN (b). The plots 
for cm(q + G , )=l/e CJ Q(q) with £gg(?) being the dielectric 
functions calculated with constrained RPA are illustrated as 
a function of \q + G\. We see that, in the limit q + G— >0, 
£m(<7 + G) converges to the two finite values. It should be 
noted here that the metallic screening effects are excluded in 
the constrained RPA calculations. In k-NCS, the higher value 
of 5.3 is the constrained-RPA dielectric constant for the elec- 
tric field perpendicular to the two-dimensional layer, while the 
smaller value of 3.8 is the in-plane value. This is character- 
istic feature of the anisotropic system. For k-CN, the higher 
dielectric constant is 4.8 and the smaller one is around 3.4. 

The data quality of the dielectric function critically af- 
fects the accuracy of the screened Coulomb interaction of W 
[Eq. ©]. In general, the dielectric function severely depends 
on the number of the conduction bands considered in the cal- 
culations B32I1 . In the ET compounds, the conduction states 
were found to be fully continuous and highly-entangled. So, 
the convergence of the dielectric function of this system was 
checked by increasing the number of the conduction bands. 
We found that the convergence requires more than 350 con- 
duction bands, which corresponds to considering the excita- 
tion up to 16 eV above the Fermi level in the polarization cal- 
culation. 
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We next show in Fig. [2] calculated ab initio macroscopic 



FIG. 2: Calculated macroscopic dielectric functions of k-(BEDT- 
TTF) 2 Cu(NCS) 2 (a) and k>(BEDT-TTF) 2 Cu 2 (CN) 3 (b) as a function 
of |q + G\. Calculated screened Coulomb interactions of k-(BEDT- 
TTF) 2 Cu(NCS) 2 (c) and k>(BEDT-TTF) 2 Cu 2 (CN) 3 (d) as a function 
of the distance between the centers of maximally localized Wannier 
orbitals. The red and green dots represent the bare and screened in- 
teractions, respectively. The solid and dotted curves denotes 1/r and 
l/(Ar), where a decay constant A was determined by the fitting to 
the ab initio data (see TABLE^. 

The lower two panels of Fig.|2]plots our calculated ab initio 
screened Coulomb interaction Vy [Eq. ([3])] with constrained 
RPA, denoted by green dots, as a function of the distance be- 
tween the centers of the MLWOs; r=\((f>i\r\(j)i) — (4>j\r\4>j)\. 
The panels (c) and (d) show the results for k-NCS and k- 
CN, respectively. The Vij decays as an isotropic function of 
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l/(Ar) (dotted line) beyond the nearest-neighbor distance (>7 
A), where A is a decay constant determined by the fitting pro- 
cedure to the ab initio data. We found that the resulting A 
is 5.1 for k-NCS and 5.0 for k-CN. For comparison, we also 
plot bare Coulomb interactions as red dots, which should de- 
cay as 1/r (solid line) 113411 . We see that the screening by the 
constrained RPA reduces the bare onsite Coulomb repulsions 
to the values less than a quarter. The onsite Hubbard U thus 
obtained is found to be 0.83 eV for k-NCS and 0.85 eV for 
k-CN. 

Our derived low-energy models differ from the conven- 
tional ones considered so far for the K-type compounds. The 
ab initio calculations give a larger ratio U / i~ 12- 15 than the 
Hiickel results ^7-8 0, & @]. Also, the ratio of the near- 
est neighbor repulsion V to the onsite Hubbard U is 0.45 
for k-NCS and 0.43 for k-CN, which are definitely larger 
than the cases of transition-metal oxides [for example, SrVC<3 
(V/U ^0-2)1, an iron-based oxypnictide LaFeAsO (^0.25) 
1 2411 . and a sodium-cluster loaded sodalite (^0.22) M26I1 . On 
top of that, we observe a remarkable long-range tail decaying 
as l/(Ar). We note that numerical studies for simple Hubbard 
models at U jt^YL without intersite repulsions are well inside 
the insulating region |4j, |6j] contrary to the metallic behavior 
of k-NCS. The results strongly suggest a relevant contribu- 
tion of the offsite interactions to the low-energy physics of the 
K-type compounds, requiring careful reconsiderations for the 
effective interactions of the K-type compounds. 

To conclude, effective low-energy Hamiltonians of organic 
compounds are established from first principles with the pa- 
rameters shown in TABLE U for k-NCS and k-CN. The de- 
rived parameters indicate that (i) the geometrical frustration 
parameter \t'/t\ is substantially smaller than the extended 
Hiickel results and k-CN estimated at |t'/t|^0.8 has turned 
out to be away from the right triangular structure and (ii) 
the onsite Coulomb repulsion (C/^0.8 eV characterized by 
U/t^l2-l5) is unexpectedly large compared to the Hiickel 
estimate given by [7/^7-8, while the intersite Coulomb in- 
teraction was found to be also appreciable. Reexaminations 
on the low-energy physics of the k-ET compounds are desired 
on the present quantitative and reliable basis, which will need 
future studies using accurate solvers for the realistic model. 
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